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Abstract 



Using the Landauer formulation of transport theory, we predict that di- 
electric quantum wires should exhibit quantized thermal conductance at low 
temperatures in a ballistic phonon regime. The quantum of thermal con- 
ductance is universal, independent of the characteristics of the material, and 
equal to ir 2 k^T/3h where ks is the Boltzmann constant, h is Planck's con- 
stant and T is the temperature. Quantized thermal conductance should be 
experimentally observable in suspended nanostructures adiabatically coupled 

to reservoirs, devices that can be realized at the present time. 
PACS: 85.30.Vw, 73.23.Ad, 63.22.+m 
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During the last two decades, the physics of electron transport in one dimension has 
attracted a great deal of attention. Some remarkable associated phenomena have been the 
quantum Hall effect discovered by von Klitzing, Dorda and Pepper |I[ and the quantized 
conductance of ballistic point contacts discovered by van Wees et al. [0 and Wharam et al. 
H Here the signature of one- dimensional conduction has been the quantization of the two- 
terminal and Hall electrical conductances in multiples of the fundamental quantum e 2 /h, as 
has been understood within the framework of Landauer theory, [f|] and of Biittiker-Landauer 
theory, || respectively. 

One-dimensional phonon transport should also be possible. However, despite the long 
standing theoretical interest in this topic that goes back to the 1920's, || the question 
whether the phonon thermal conductance should be quantized in one dimension has to 
our knowledge not been addressed either theoretically or experimentally. Recent advances 
in nanotechnology have made experimental investigation of this question feasible; we note 
especially the detection of nanowire phonon subbands by Seyler and Wybourne and the 
measurement of the thermal conductance of a suspended nanostructure by Tighe, Worlock 
and Roukes [§]. 

The purpose of this Letter is to demonstrate theoretically that in a low temperature 
regime dominated by ballistic massless phonon modes the phonon thermal conductance 
of a one dimensional quantum wire is quantized, the fundamental quantum of thermal 
conductance being ir 2 k^T/3h, where /eg is the Boltzmann constant, h is Planck's constant 
and T is the temperature. We also establish the conditions that should be met for the 
experimental observation of this novel phenomenon. 

Our starting point is the Landauer energy flux 

f°° dk 

Q = J2 7T~ tWa(k) V a (k) (r) R - T) L ) C a (k) (1) 

carried by a quantum wire connecting two reservoirs labeled R and L. Here u a (k) and v a (k) 
are the frequency and velocity of normal mode a of the quantum wire with wave-vector k, 
( a (k) is the phonon transmission probability through the wire and rji{uj) = 1/ {e 1Vj} ' kBTi — 1) 
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represents the thermal distribution of phonons in the reservoirs, assumed to be a Planck 
distribution at temperature Tj. The cross sectional area of the wire is assumed to be of 
the order of hundreds of nm 2 , so that the lateral confinement produces finite gaps in the 
dispersion relation of the phonon frequencies. Eq. (TfJ) transforms to 

1 f°° 

Q = —^l du hw (t) R (uj) - r) L (uj)) ( a (uj) (2) 

since the phonon velocity v a (k) = dui a /dk is canceled by the ID density of states g(u a ) = 
dk/duj a . 

The reservoir-to-reservoir thermal conductance of the wire is k = Q/AT, where AT = 
Tr — Tl is the temperature difference between the reservoirs. It follows that 

We have separated Eq.(H) in two parts, so that the first term represents the conductance 
of the massless modes with uj a (0) = and the second is the contribution to the thermal 
conductance due to the higher energy modes, which have a finite cut-off frequency ov(0) ^ 0. 

For the moment let us assume perfectly adiabatic contact between the thermal reservoirs 
and the ballistic quantum wire, so that £(a>) = 1. For this idealized case integration of 
Eq.(|) yields 

kin 2 /T R + T L \ , r kl in 2 /T R + T L \ 1 r - , , /T , - /T . n ] 

where x^lk = (0) / ksT^L) , f° r the modes with finite cut-off frequencies. 

Equation (£|) predicts a remarkable behavior for the thermal transport properties of ID 
phonon systems: Each massless mode presents a universal thermal conductance equal to the 
product of the averaged reservoir temperature and the universal constant k 2 B n 2 /3h. The 
higher energy modes, on the other hand, show a dependence on the intrinsic properties of 
the material and on the geometrical parameters of the sample through the cut-off frequencies 
uv(0). However their contribution to k is exponentially small at low temperatures. In the 
limit AT — > we find 



p 2 L2 ^a' f_2 2 x -\ 

K = ^TN a + %T V - + /(z ) + (5) 
3/i h \ 3 17 V ; e x ° - 1 J V ; 

with /(x) = 2dilog(e x ) and xo = hu a i(0)/kBT. 

Expressions (^) and (|5|) represent an idealized case in which the transmission of phonons 
from one reservoir to the other happens without reflections. A more realistic model has 
to incorporate the effects of the reflections caused by the contacts between the reservoirs 
and the ID wire. The remainder of this paper is devoted to the study of realistic physical 
systems in which the universal thermal conductance of phonons could be observed. The 
systems we consider are similar to the experimental device of Tighe, Worlock and Roukes 
|| where a quasi-lD quantum wire connects quasi-2D reservoirs. 

We begin our analysis by examining the dispersion relations of low energy phonons in a 
quasi-2D system. At low temperatures of the order of a few degrees Kelvin, the dominant 
phonon wavelength is much larger than the lattice parameter and the model of an elastic 
continuum can be used. We assume an isotropic crystal, which can represent Si or GaAs 
at low temperatures to a good approximation. The optical phonons are not considered 
because of their high energy. In this model Rayleigh-Lamb modes || describe the quasi-2D 
acoustic phonons. Figure [3] shows the frequency spectra of the lowest energy symmetric (solid 
lines) and anti-symmetric (dashed lines) modes of a quasi-2D system of 50nm thickness. 
The parameters used are those of GaAs. The symmetries of the modes refer to atomic 
displacements in the z direction, perpendicular to the quasi-2D system. The energy scale 
is degrees Kelvin, so that one may estimate the temperature at which the modes with 
finite cut-off frequency begin to contribute to thermal transport. These cut-off frequencies 
increase as the sample thickness is reduced. Together with the 2D modes of the reservoirs 
we have plotted in Fig. |I| the lowest energy branches of the longitudinal (circles), transverse 
(squares), shear (diamonds) and torsional modes (triangles) of a quasi-lD quantum wire 
with a square cross section of 50nm x 50nm. The fundamental torsional mode is given by 
the simple beam theory, which is well justified by the fact that numerical calculations [ 10[ 



show that its phase velocity is practically constant for all frequencies. The transverse and 



shear modes are coupled by the Timoshenko equation JTTJ. These quasi-lD modes agree very- 
well at low energies with the 2D modes of the reservoirs. Therefore they should propagate 
readily from the reservoirs to the wire. Since expressions (f|) and (|5|) do not depend on the 
phonon velocity, but only on the cut-off frequencies, the modes shown in Fig. [1| represent 
the system very well up to temperatures of the order of IK. 

Nevertheless, there still remains the question of how reflections would change the results 
given in Eq. ([|) and (f|). We will discuss this here in detail for the case of the longitudinal 
modes. For a quasi- ID wire whose cross sectional dimensions are much smaller than its 
length, the transmitted modes are well described by the normal modes of a long beam. Let 
us suppose that the cross section of the beam varies along its length, assumed to be in 
the x direction. The equation of motion of a longitudinal plane wave traveling along the x 
direction is given by [|ll|] 

d 2 u 1 dA(x) du 1 d 2 u 
dx 2 A(x) dx dx vf dt 2 

A(x) is the cross sectional area of the beam and V\ is the velocity of the longitudinal mode 



vi = yY/p, Y is Young's modulus and p the density. To solve Eq.(|]) we must specify 
A(x), so that the transmission coefficients ((w) are determined by the shape of the contacts 
between the reservoirs and the ID wire. We consider two contact geometries: conical, for 
which A(x) = A Q tan(6)(x + x'), and catenoidal where A(x) = A cosh 2 (x/X). These are 
illustrated in Fig. |2| by schemes (A) and (B), respectively. 9 is the angle of flare of the 
conical contact and A is the characteristic length of the catenoid. The broad regions at 
both extremes of the structures are 40 times wider than the ID channel (50 nm wide). The 
thickness (in the direction perpendicular to the plane of the figure, z) is constant along the 
whole structure and chosen to be 50nm. Considering initially the case of conical contacts, 
two distinct wave equations result for the straight (I) and conical (II) parts of the structure 

d 2 ui 1 d 2 ui 



dx 2 vf dt 2 



(7) 



d 2 u n | 1 du H _ 1 d 2 u n 
dx 2 x + x' dx vf dt 2 
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with solutions 

ui{x,t) =Be i{±kx ~ wt) (9) 
unix, t)=C [J Q (k(x + x')) ± iN (k(x + x'))} e~ iwt . (10) 

These solutions are matched at the interfaces between the conical and straight regions, 
where the continuity of the particle velocities u(x,t) and stresses a = A(x)Y 9u q^ are 
required. These boundary conditions lead us to a system of coupled equations that is solved 
for the transmission coefficients of plane waves traveling along the structure. FigfJ shows the 
transmission coefficient ((h) as a function of the longitudinal wave- vector k of a structure 
with conical contacts in which 9 = tt/6 (solid line). The peaks are resonances associated with 
reflections at the ends of the conical regions and of the ID wire. For k = the transmission 
coefficient is equal to one. The overall behavior of ((k) is determined by the conical shape 
of the contact, that guarantees a finite transmission for all frequencies, but is inefficient for 
a wide range of frequencies. 

The catenoidal contacts present a better response for lower values of k. In this case the 
equations that define the transmitted wave are 

d 2 ui 1 d 2 ui 



dx 2 v? dt 2 



(11) 



d 2 un 2tanh(x/X) dun 1 d 2 u n , . 



dx 2 A dx vf dt 2 



with solutions 



Ul (x,t) = Be i{±kx ~ wt) (13) 
u u (x,t) = Ce i{±kx - wt \cosh(x/\))- 1 . (14) 

The transmission coefficients obtained in this case are illustrated by the dashed curve in 
Fig. 0. Substitution of the function ujj in Eq. ( p~2|) results in the frequency spectrum oj 2 = 
vf(k 2 + (1/A) 2 ) for a wave traveling along the catenoidal contact, which has a cut-off at 
frequency u = vi/X. For u < u the wave is evanescent, but a resonance guarantees unitary 
transmission for k=0. However the cut-off frequency becomes smaller as the parameter 
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A that characterizes the length of the catenoid increases. The efficiency of transmission 
of longitudinal waves in this system is also limited by reflections that happen on passing 
between the straight region and the curved one. These reflections can be avoided for a 
catenoidal contact that is not laterally limited. This type of contact is exemplified by the 
structure (C) in Fig. |2|, the infinite catenoidal contact. The transmission ((k) for this case 
is 1 for all wave vectors k, but the cut-off at ujq = vi/X still affects the overall transmission. 

With a knowledge of ({uj) it is possible to obtain realistic results for the ID thermal 
conductance. Initially we analyze the case of ideal transmission for all modes, which is a 
very good approximation for the large k acoustic waves that dominate the thermal transport 
between lOOmK and IK. Six distinct modes contribute to the thermal transport through 
the wire below IK: a longitudinal mode, 2 transverse modes, a torsional mode and 2 shear 
modes. Among these only the shear modes present a cut-off (huj Q = 1.62K) that is the result 
of lateral confinement (see Fig. p]). It was verified that the higher branches of these modes 
do not contribute significantly to the thermal transport at temperatures lower than IK. We 
considered the limit AT — > of Eq.(|5|) and calculated the ideal thermal conductance divided 
by temperature of the quantum wire. The result is shown in Fig. |3| (refer to the left scale). 
For the ideal contact, at low temperatures, this yields the quantized thermal conductance 
k%ir 2 /3h times the number of massless modes (N a =4). As the temperature approaches IK 
the shear modes begin to contribute to the process (increasing the thermal conductance) 
and the plateau that is the signature of quantized thermal conductance terminates. Now 
considering the transmission coefficient of the longitudinal mode, we have plotted its con- 
tribution to the thermal conductance for various contact shapes (refer to the right scale of 
Fig. |3|). In this case the plateau is modified at the low temperature side by reflections at the 
contacts and, to a smaller extent, by the cut-off of the frequency spectrum of the catenoidal 
contacts. Nonetheless for the catenoidal contacts a very distinctive plateau characteristic of 
quantized thermal conductance is clearly visible over a wide range of temperatures. There- 
fore based on the behavior of the longitudinal mode, a signature of the quantized thermal 
conductance of dielectric quantum wires can be expected between 30mK and 300mK. We 
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note that the principle of adiabatic matching that underlies our discussion of the longi- 
tudinal mode is very general and applies to the other massless modes as well. Therefore 
these should exhibit qualitatively similar behavior, including universal quantized thermal 
conductance over a similar range of temperatures. Since adiabatic wave propagation does 
not require an isotropic medium, our predictions apply even if there is phonon focusing due 
to crystal anisotropy. Similar results also apply to wires with cylindrical symmetry which 
admit analytic solutions for all of the modes, but will be more difficult to fabricate. 

In this study we did not consider the effects of surface roughness f| and other defects 
that may limit the transmission of phonons through the quantum wire. However these 
are technological limitations that should be overcome with improved control of the growth 
process. They should also be of less concern at lower temperatures. |12| Furthermore, studies 
of ID electron transport have shown that quantization plateaus can exist in the presence of 



defect scattering under certain conditions. |L3|. 

In summary, we predict the existence in ID systems of a universal quantum of thermal 
conductance due to phonons that is equal to k B ir 2 /3h times the temperature, for the lowest 
energy modes. The conductance of the higher energy modes is influenced by the geomet- 
rical and intrinsic parameters of the system through the gaps in the frequency spectrum. 
Realistic structures were investigated and it was found that contacts with catenoidal shape 
should be those that better represent the ideal case in experiments. Finally, we note that 
the Wiedemann- Franz law applies in ID fll4| , [|15f and has been tested experimentally for 
ballistic point contacts. [|16j It predicts |15| that the ID thermal conductance of ballistic 
electrons should also be quantized in multiples of k 2 B ir 2 T '/3h. Thus we arrive at the surpris- 
ing conclusion that the low temperature ID thermal conductance of ballistic phonons and 
electrons is described by the same universal quantum even though bosons and fermions obey 
different statistics. 

We thank Michael Roukes for stimulating discussions. This work was supported by 
NSERC of Canada. 
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FIGURES 

FIG. 1. Lowest energy acoustic modes of a quasi-2D system of thickness 50nm and of a long 
wire of square cross section 50nm x 50nm. Thick solid (dashed) lines represent the symmetric 
(anti-symmetric) modes of the quasi-2D system. Thin lines with circles, squares, diamonds and 
triangles represent the longitudinal, transverse, shear and torsional modes of the wire, respectively. 
The elastic parameters are for GaAs. 

FIG. 2. Transmission coefficients vs. longitudinal wave- vector k. A) conical contact, 6 = ir/6; 
the quasi-lD wire at the center of the structure has cross-section 50nm x 50nm and length 1//. 
B) finite catenoidal contact of characteristic length A = 0.86//. C) infinite catenoidal contact of 
A = 0.86//. The straight portions at the ends of structures (A) and (B) have width 2//. 

FIG. 3. Left scale: Thermal conductance of a quantum wire with ideal contacts divided 
by temperature. Right scale: contribution to thermal conductance due to the longitudinal mode 
for various contact shapes: infinite catenoid for A = 4.6// (solid), finite catenoid for A = 4.6// 
(dot-dashed), catenoid with A = 0.86// (long-dashed) and conic for 6 = ir/6 (solid with circles). 

< d eal = kWT/Zh. 
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